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Discovery of shocked H2 around OH 231.8+4.2 
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ABSTRACT 

Wc present X-band integral field observations of the cireumstellar envelope of 
the evolved star OH 231.8+4.2. Spatial and spectral information were simultaneously 
acquired using the SINFONI integral field unit, with adaptive optics, on the Very Large 
Telescope. The observations reveal the discovery of H2 emission (1) around the centre 
of the nebula and (2) located in clumps along the Western side of the Northern lobe, 
presumably associated with the strong shocks that stimulate the previously reported 
Ha emission at the same location. An observed H2 1— >0/2— !'1S(1) line ratio of 8.3±1.9 
was calculated for the central field, a value consistent with shock excitation. 

Key words: molecular data - cireumstellar matter - stars: AGB and post-AGB - 
stars: individual(OH 231.8+4.2) - shock waves 



1 INTRODUCTION 

OH 231.8+4.2 (hereafter OH231) is an O-rich late spectral 
type (M) central star (Mira variable, QX Pup) with bipo - 
lar high- velocity dust and gas outflows (jAlcolea et al.|[200lh , 
filamentary structures observed in scattered and molecular 
line emission, and large angular size (10"x60"). Often la- 
beled as a post-AGB object or pre-planetary nebula, the 
presence of both a Mira central s tar and a main-sequence 
comp anion of spectral type A (jSanchez Contreras et all 
l2004h seems to contradict this classification. OH231 is more 
likely a D-type bipolar symbiotic system (|Frew fc Pajkej 
I2OIOI ). However in some cases, morphological similarities do 
exist between post-AGB and symbiotic objects, most strik- 
ingly the presence of highly coUimated bipolar nebulae. It is 
via fast-coUimated outflows that these stars shape their sur- 
rounding nebula. Understanding the development and origin 
of these fast outflows is critical for advancing hyd rodynam- 
ical m odels of wind interaction. Recent work by iLee et al] 

attempting to reproduce the high velocity molecular 
emission in AFGL 618 using coUimated fast wind models, 
emphasises the need for further observations and model de- 
velopment in this area. 

OH231 has been the subject of many s tudies span- 
ning m ultiple wavelength ranges, for example: ICohen et al.l 

were first to propo se the existence of a binary com- 
panion; iBuiarrabal et al" 

I (I2OO2') imaged the shape of the 
shocks using Ha (reproduced in Fig. [T] [a]) detected w ith 
the Huhble Space Telescope (HST); iMeakin et all (j2003l ) re- 
port HST/mCMOS NIR images of the dust distribution and 
hence a high resolution map of the extinction through the 
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nebula. iMatsuura et all (|2006l ) using the MIDI and NACO 
instruments on the Very Large Telescope (VLT) detected 
a compact cireumstellar disc. The envelope of OH231 is 
also known to be rich in molecular species (e.g. H2O, OH, 
and SiO) however previous studies in the NIR have all re- 
turned null detection s of H2 (e.g. IWeintraub et al] Il998l : 
iMatsuura et al.ll2006t ). 

In this Letter, we present the results of preliminary ob- 
servations of OH231 at NIR (if-band) wavelengths show- 
ing for the first time the presence of shock-excited H2. 
Throughout this work we assume OH231 is a member of the 
open cluster M46 at a distance of 1.3 kpc l|jura fc Morris! 
Il985h . The origin of the coordinate system used in Fig- 
ures 1,3, and 4 is given by the location of the SiO 
maser e mission at RA=07''42"'1 6^93, Dec=-14°42'50".2 
(J2000) (jSanchez Contreras et all[2002 ). and the inclination 
angle of the bipolar axis is 36° to the plane of the sky 
l|Kastner et al.lll992l ). 



2 OBSERVATIONS AND DATA REDUCTION 



The data were taken using the SINFONI (lEisenhauer et aP 
I2OO3I ) instrument located on UT4 at VLT at Paranal, Chile, 
on the l^VS'"* Feb 2010. We use the lowest resolution mode 
(LRM), corresponding to the widest field-of-view (8" x 8") 
with adaptive optics (AO) and a natural guide star (NGS). 
A plate scale of 250x125 mas pixel"^, and a spectral and 
velocity resolution of 4580 and 66 km s~^, respectively, are 
available at this resolution (for a dispersion of 2.45 A/pix 
and line FWHM of 1.96 pixels). All observations utilised 
the Tf-band (2.2 ^m) filter which covers many ro- vibrational 
H2 emission lines. The ambient seeing varied from ~ 0.6" 
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Table 1. Details of the VLT/SINFONI observations for OH231 and 
the telluric standard stars taken on f^*/ 2"*^ February 2010. 





OH231 




HD 


South 


Central 


North 


75004 63487 


Tot. Exp. (sees) 600 


600480 


800 


1 10 


K (mag) 






7.268t 7.674t 


Night 1 


1 2 


1 


1 2 



t 2mass magnitude from; SIMBAD 



to 1.1" during the observations. The OH231 observations 
consisted of three fields focused on (1) the edge of the 
Northern lobe, (2) the central region, and (3) the middle 
of the Southern lobe (labelled N, C, S in Fig. [1] [b]). No 
H2 was detected in the Southern field and will not be dis- 
cussed further. Table [l] summarises exposures for each of the 
three fields. Telluric standard stars used for calibration are 
HD 75004 (GOV), and HD 63487 (G2V) for night one and 
two, respectively. An average AO-corrected PSF of ~340 
mas FWHM is estimated from the standard stars. The data 
were reduced using the ESO common pipeline library to a 
wavelength-calibrated datacube and further analysed using 
both PyRApQ. Wavelength calibration errors were corrected 
by comparison of OH emission lines with a high resolution 
template. Quoted velocities were adjusted to local standard 
of rest (LSR) velocities using 22.88 and 23.30 km s~^ cor- 
rections for night one and night two, respectively. Line maps 
(Fig. [1] [c,d]) were generated by fitting the H2 emission lines 
with a Gaussian profile; initial fit parameters (FWHMq, cen- 
tral wavelengtho, etc.) were determined by manually fitting 
an individual H2 fine, only fines with a FWHM « FWHMq 
make up the final line maps. Signal-to-Noise (S/N) of the 
line maps was enhanced by smoothing the data with a 
2x2 pi xel boxcar filter. Line rest w avelength information 
is from lBlack fc van Dishoeckl (|l987l ). 



3 RESULTS 

We report the detection of several H2 emission lines arising 
from both the centre and Northern lobe of OH231. In Fig- 
ure [1] (c,d) we present continuum-subtracted line maps of 
the 1— >-OS(l) transition for both fields, showing clearly H2 
emission arising from the centre of OH231 (Fig. [T] [d]) and 
from knots of material in the Northern lobe (Fig. [U [c]). In 
Fig. [T] (a) the contours of the detected H2 are shown in re- 
lation to Ha emission, while Fig. [1] (b) gives the location of 
the H2 relative to the strong continuum emission as imaged 
with a 2.12 ^im filter. An integrated spectrum of the Western 
region of the centrally located H2 is shown in Figure (2] from 
which we note, a) several S- and Q-branch ro- vibrational H2 
lines, and b) a CO bandhead absorption feature (~2.3 \im). 

Channel maps, extracted from the SINFONI datacube, 
are presented in Figures[3]and|4]for the central and Northern 
regions respectively, showing how the distribution of the H2 
changes across the line profile. 

Table [2] presents the fiux measurements for the central 
and Northern fields, for lines with flux errors less than ~50 

^ PyRAF is a produet of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 



Table 2. Line fluxes, F, (in units of IQ-^^ W m-^) of the ob- 
served H2 lines for the central and Northern fields. Measured peak 
line wavelengths (in (xm) are given for the central field. 



Line 


Central 


Northern 




F 


F 


1~!>0S(3) 


1.9579 


67.2 ± 14.2 


11.5 ± 1.10 


1^0S(2) 


2.0342 


9.71 ± 3.36* 




l^OS(l) 


2.1222 


91.5 ± 6.16 


9.69 ± 0.76 


2-^lS(l) 


2.2481 


7.73 ± 3.38* 




l^OQ(l) 


2.4071 


46.4 ± 14.5 


11.9 ± 2.95 


1^0Q(3) 


2.4241 


42.2 ± 24.0 


13.3 ± 4.42 



* Flux measured for Western H2 region only. 

* Measurement confined to elliptical region (Fig. [T] [d] ) . 



per cent. Accurate 2^^1S(3), 1-5>0S(0), and 1^.0Q(2) line 
flux measurements are not possible due to the presence of 
strong sky subtraction residuals at these wavelengths. We 
confine the calculation of the 1— >-0/2— >-lS(l) ratio for the 
central field to the area marked by the ellipse in Fig. [T] 
(d), where the 2^1S(1) flux error is smallest. In this re- 
gion we calculate a 1—^0/2— >-lS(l) ratio of 8.3±1.9 prior to 
extinction correction. We do not detect any 2— >-lS(l) flux 
in the Eastern H2 region, instead we estimate an upper 
2-^lS(l) limit (3cr) for the fiux in this region of 1.3x10"" 
W m~^, which in turn places a lower limit of ^^9.5 on the 
1^-0/2— >-lS(l) ratio for this region. 



4 ANALYSIS AND DISCUSSION 

The 1— >-0/2— >-lS(l) ratio is a useful discriminator of ex- 
citation mechanisms. Pure fiuorescence will yield a value 
of ~ 2 while a value ~ 10 indicates the excitation of the 
gas is being driven by shocks. However it is worth not- 
ing, these values depend on shock velocity and pre-shock 
gas density, with shocked H2 capable of producing values 
as low as «4 (|Smithlll995l ). while fluo rescence can produce 
value s approaching those of shocks (|Hollenbach fc Nattal 
'1995"). The 1^0/2^18(1) ratio values, given above, sug- 
gest that shocks might be the main excitation mechanism, 
which agrees with the detection of shock-excited HCO^ in 
t he centre of OH231 noted i n the position- velocity diagrams 
of lSanchez Contreras et al] (|2000l V In order to determine the 
intrinsic 1— >0/2-^-lS(l) ratio, it is necessary to remove the 
effects of extinction. In the K-hanA, it is sometimes possi- 
ble to derive the level of extinctio n via the comparis on of 
S-/Q-branch H2 emission lines (see lSmith et al.ll2003l ). Un- 
fortunately, due to poor atmospheric transmission above 2.4 
\im, it was not possible to derive a sensible estimate for ex- 
tinction using the 1— >0 S(l) and Q(3) lines. 

Extinction values of Ak=3-4 (mag) for the central 
region are e stimated from the {K — L') colour map of 
OH231 (Kast ner et"al]|l998l '). Although, most likely an over- 
estimate, adjusting the Western 1— >-0/2— >1S(1) ratio for 
these levels of extinction yields, for example, a value of 
11.7±2.6 (j4k=4). It is clear that any adjustment for extinc- 
tion will increase the observed 1^-0/2— >-lS(l) ratio, pushing 
it further towards the shock regime. 

Using a typical value of 5.3x10"^^ niag cm^ for the 
extinction per unit column density of hydrogen, Ay/NiYi.), 
one can estimate the hydrogen column density implied by 
an Ak=4; yielding a iV(H) for the H2 regions of ~7.2xl0^^ 
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Figure 1. Narrow- band images and line maps for OH231; a star symbol indicates the SiO maser position (see text), inset labels display 
the telescope/instrument information, (a) Continuum-subtracted Ha image with the 1— >'0S(1) H2 contours overlaid in white, (b) 2.12 i^m 
image, showing mainly scattered light, with H2 contours overlaid, SINFONI North, Central, and South (labeled N, C, S) fields are marked 
by green squares, inset compass shows orientation, (c) 1— >0S(1) H2 line map showing the extent of the emission located along the edge of 
the Northern lobe, two knots (A, B) of H2 emission are indicated by arrows; superposed contours indicate the location of the Ha emission, 
(d) 1— ^OS(l) H2 line map for the central region showing the position/extent of the region used for 1— >'0/2— >'1S(1) calculation (ellipse). 
Colour bar units are in W m~^ px~^. T he Ha data, from the H ubble Legacy Archive (see Acknowledgm ents), has been c ontin uum- 
subtracted using the procedure outlined in lBuiarrabal et al. I1 I2OO2I) . The 2.12 \im data was first published in lMatsuura et al] 1I2OO6I ) and 
has been reprocessed using the ESO pipeline. 



cm ^ . Using a typical column length of the H2 emitting 
region s of l"(~2.0xl0^® cm), as used in iBuiarrabal et al.l 
we estimate an average density for the H2 regions 
of 3.5x10^ cm~'^. This value is in good agreement with 
lAlcolea et al.l (|200ll ) who estimate an average central density 
of ~3.0xl0® cm~^. The H2 emission is most likely originat- 
ing in the dense equatorial regions surrounding the central 
star. 

4.1 Equatorial Region 

The l^'OS(l) H2 hue map (see Fig. [T] [d]) shows two re- 
gions of H2 oriented at a position angle (PA) of 113°, from 
brightest to faintest peak. If we assume that the distribution 
of the H2 around the central region of OH231 is in a disc 
configuration then by measuring the major and minor axis 
we can estimate the inclination angle of the equatorial disc 
with respect to the plane of the sky from the H2 data. Mea- 
surements for -Rmax/iimin are determined by superposing a 
full ellipse onto the H2 line map. We find -Rmax= 4.06" and 
J?min= 2.76"from the centre of the ellipse, and using the re- 
lation i = sin~^(_Rmin/i?max), we find i =43°±8. This result 
is in good agreement with previously published values. 



The location and orientation of this H2 adds to the al- 
ready complex picture of the equatorial region of OH231. 
Some previously reported structures, from smallest to 
largest, include: 

- an equatorial torus of SiO maser emission (~2R*) that 
migh t lie on the innermost edge of an expanding SO 
disc jSanchez Contreras et"al] [2000. 2002); 

- a centrally located compact disc of circumstcUa r mate- 
rial w ith inner R=0.03-0.04"/40-50 AU (Matsuur aet all 
iooi); 

- an opaque flared disc with outer R=0.2 5"/330 AU re- 
vealed in mid-IR images jjura et al.ll2002l ); 

- a slowly expanding disc with characteristic 
R=0.9"/116 AU, detected via the S O (J=22->li) 
transition (Sanchez C ontreras et al.ll2000l ); 

- a torus of OH maser emission with outer R=2.5"/3250 
AU (fZiilstra et al. 2001); 

- an expanding hollow cylinder of HCO^ with a 
characteristic radius equal t o the OH torus radius 
jSanchez Contreras et aLll200(j ). 

- a halo of scattered light at 

R~4"/5200 AU jMeakin et al.ll2003l ): 

From our observations the geometry of the H2 region is un- 
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Figure 2. An integrated spectrum of the Western side of the 
central H2 showing the detected hnes, the position of the Br7 
recombination hne is shown for reference. Extracted spectra lo- 
cation marked in Fig. [l](d). 



clear, however we offer two possibilities: (f) A disc of H2: 
Figure [1] (d) shows what might be interpreted as an in- 
complete disc of H2, which fits with the series of concen- 
tric disc/tori structures listed above. To understand why 
we observe only emission from the edges and not the front 
of the disc, a comparison of the noise in the continuum 
at the front of the putative disc, to the amplitude of the 
1— >-OS(l) line peak shows both to be of the same order. 
We might then attribute the 'missing' £[2 in this region to 
variations in the continuum. We would not expect to ob- 
serve the back of the disc due to the high level of extinction 
through the nebula. (2) A shell of H2: if the H2 is situated 
in an axisymmetric shell surrounding the central star, we 
might explain the geometry of the H2 regions by assuming a 
density contrast between the poles and equator. This, com- 
bined with an increased column depth at the edge of the 
shell would manifest itself as t wo arcs of H2 emission sit- 
uated equatorially (as noted by iLowe fc Gledhil]|[2005l . for 
IRAS 19306-1-1407). This is supported by the detection of 
a shell of higher density gas a nd dust at the same location 
as the H2 (|Meakin et al.ll2003l ). Both scenarios offer plausi- 
ble explanations for the geometry of the H2 emitting region, 
however further observations are needed in order to favour 
one. 

We fit the 1->0S(1) line profile yielding a Vlsr=36±17 
km s^^and FWHM^IOO km s^^for both Eastern and West- 
ern regions of H2 emission, a value consistent with the sys- 
temic velocity. Channel maps of the l^OS(l) line (Fig. ^ 
show no significant change in the distribution of the H2. 

In an attempt to explain the lack of re ported H2 in this 
object , we note two previous studies: (1) IWeintraub et al] 
l|l998l) give a Sa upper limit of 10~^ ergs cm"'^ s~^ ster~^ 
(Icr limit =6.6x10"^* W m"'^) for the surface brightness 
of the 1— >-OS(l) line towards OH231. This places their mea- 
surement limit close to the 1— >0S(1) line strength given in 
Table [2j implying that the 1— >-0S(l) line in their observa- 
tions would be difficult to distingush from noise, or pos- 
sibly that their slit position, aligned East- West across in- 
ferred central star pos ition, did not include the H2 sites; (2) 
iMatsuura et al.l(|2006l '). using the ISAAC instrument on VLT, 
do not report any detection of H2. However, this can be ex- 
plained due to the orientation of the slit, aligned from South- 
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Figure 3. 'Channel maps' of the four central channels of the 
1— ^OS(l) line from the central region of OH231. Inset labels 
show datacube channel and corresponding Vlsr velocity. Note: 
in ch794, the central region (marked by a green box) is residual 
from the continuum subtraction and not true H2 emission. 



West to North-East along the major axis (private comm.), 
with a slit- width of 0.8", i.e., the central £[2 emission site 
was not covered. 



4.2 Northern Region 

Figure[T](c) shows the line map for the 1— >0S(1) transition 
in the Northern region. Most notable are the two knots of 
H2, labeled A and B for the top and bottom knot, respec- 
tively. The Ha emission contours are superposed on the H2 
line map. The peak intensity of the H2 emission knots are 
slightly offset from the two Ha emission knots (Fig. [U [c]), 
however this small offset can be accounted for by the motion 
of th e outflow (e.g.l/~150 km s~^, from [Suiarrabal et al.l 
2002) during the 10 years between both sets of observa- 
tions. It is most likely that the optical and NIR line emission 
arise from the same shock event. Weaker 1— >-OS(l) emission 
is noted tracing the Ha edge of the bipolar outflow in the 
lower portion of the 1->0S(1) line map. The NACO 2.12 \im 
image does not show any trace of H2 emission in this region 
(see contours in Fig. [U [b]). 

The l-^OS(l) line is spread over ~six spectral pixels, 
the four central channels are presented in Fig. U Examina- 
tion of the line proflle and channel maps allows us to probe 
the kinematics of the H2 in this region, revealing two main 
features: 1) knot B which is persistent in all channel maps 
2) the reduction of peak H2 intensity in knot A at velocities 
^ -56 km s~^. Both of these H2 structure s lie in the diffuse 
exten ded region (labeled Bi in Fig. 4 of iBuiarrabal et all 
l2002h perpendicular to the axial flow, with a quoted Ha ve- 
locity of 150 km s~^ which is in good agreement with not 
only our deprojected H2 v elocities (Vk^^HO km s~^) bu t 
also with HCO'*' velocities jSanchez Contreras et al1l2000l ). 
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Figure 4. 1— >-OS(l) line 'channel maps' of the four central chan- 
nels extracted from the Northern region of OH231. Knot A and 
B are regions of most intense H2 emission. Labels as in Fig[3] 



In the case of knot A, the spectral line is strongly peaked 
in a single channel, while in knot B the spectral line peak is 
spread over two channels. Fitting the line profiles of knots 
A and B, yields LSR velocities of Vlsr=-8 km s^^, and -30 
km s~^, with deprojected velocities of -75 km s^^and -110 
km s~^, respectively. This might indicate that further from 
the Ha bow apex (Fig. [T][a]), we are starting to see emission 
originating from the front and back of knot B, while knot A 
displays a narrower range of velocities, i.e., a single peak in 
its spectral profile. It is worth noting that due to the slit- 
length of ISAAC (120") coupled with the null H2 detection 
discussed in i]4.1l the II2 in the Northern region might be 
confined to the wings of the bow shock. 

The existence of fast moving shocked H2 has pre v iously 
been noted in other objects, for example, I Cox et akl (|2003l ) 
detect high velocity H2 ~220-340 km s~^ (dependent on 
adopted inclination angle) in the outflows of AFGL 618. It 
is currently unknown exactly how shocked H2 can be trav- 
elling at this speed without complete dissociation. Further 
high resolution mapping of the H2 is necessary in order to 
resolve the shock surfaces, more accurately measure the H2 
kinematics, and apply shock models to this region. 



5 CONCLUSIONS 

We have presented VLT/sinfoni integral fleld observations 
of OH231, revealing the presence of several ro-vibrational 
H2 lines. The main conclusions are: 

- The discovery of H2 emission near the centre of OH231, 
possibly located at the edge of an axisymmetric shell or 
an incomplete disc. 

- A 1^0/2-i>lS(l) value of 8.3±1.9 is found for the equa- 
torial H2, suggesting a collisional excitation mechanism. 



- Our observations of the central shell/disc of H2 show no 
velocity structure. However, higher S/N and/or velocity 
resolution data are needed to accurately probe the kine- 
matics in this region. 

- We detect fast-moving H2 (~110 km s~^, along the 
bipolar axis) via the 1— >-0S(l) transition along the North- 
western tip of the nebula, a region where a strong Hq bow 
shock exists. Due to the small FOV of our observations, 
the full extent of this H2 is unknown. 
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